function Ry = Roty(theta)
    Ry = zeros(3,3);
    Ry = [cos(theta)    0     sin(theta);
           0            1              0;
          -sin(theta)    0    cos(theta) ;];
end